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SUMMARY 


This paper presents a method of evaluation of Wiener-Hopf technique solutions for sound 
radiation from an unflanged circular duct with infinitely thin duct wall, including mean flows. 


1. INTRODUCTION 


Sound radiation from circular ducts is a classical acoustics problem. Exact solutions were 
previously reported: the Wiener-Hopf technique was used for radiation of propagating modes from 
a circular duct with negligibly thin duct wall (refs. 1-3), and the hyperboloidal wave function was 
defined and employed for radiation from duct with various types of termination, including a plane 
flange and horns (ref. 4). Exact solutions undoubtedly help one to gain physical insight into the 
problem and can often be used in practical designs. In this electronic computation age, another 
significant role of exact solutions is defined as means of cross examination of results of numerical 
techniques. These techniques, embraced as computational aeroacoustics, are just starting to attract 
widespread attention as a potential tool in attacking important aeroacoustic problems for which 
quantitative solutions are not available. 

Despite the elegance of the closed form solutions with the Wiener-Hopf technique, numerical 
presentations have been limited to mere demonstrations of its capability. As a matter of fact, no 
computer program is publicly available for its numerical evaluation. Numerical evaluation of the 
Wiener-Hopf solution is not straightforward; it requires the exercise of extreme care and, often, 
sophisticated mathematical tricks. This paper attempts to provide a comprehensive mathematical 
procedure for evaluation of Wiener-Hopf solutions. 

In section 2, acoustic waves will be briefly reviewed for in-duct propagation and radiation. In 
section 3, the Wiener-Hopf technique is applied to obtain solutions, and section 4 is devoted to the 
evaluation of integrals involved in the solutions. 


* Partly based on two consulting reports submitted to Pratt and Whitney Aircraft, May 25, 1976, and December 17, 1976. 
1 Massachusetts Institute of Technology, Cambridge, Massachusetts. 


2. REVIEW OF DUCT ACOUSTICS 


Duct acoustics will be briefly reviewed here for aspects relevant to the present problem. This 
review is also intended for clarification of terminology and nomenclature used in this paper. 

The wave equation for the acoustic pressure, p, in flow is 
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Here c is speed of sound, and V the mean flow velocity which is assumed to have only the axial 
component. The analysis is confined to a steady wave with the harmonic time dependence e lCOt , 
and axial angle dependence e tm< ^ , where m is an integer called the circumferential mode number. 
Equation ( 1) is then written for circular cylindrical coordinates (r, <j>, x ) as 
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where k = co/c and M = V/c. 

The sound radiation from an unflanged circular duct is schematically displayed in figure 1. With 
reference to this figure, the entire region is divided into two: region 1 for r < a, and region 2 for 
r> a, where a is the duct radius. The subscripts 1 and 2 will be used from now on to indicate, 
respectively, regions 1 and 2, unless specified otherwise. The mean flow velocities are assumed to 
be uniform in each region and are denoted by Vj and . For V) & Vj, there will be the mean flow 
mismatch at r = a, for x > 0. The sound speed can differ for the two regions for reasons such as 
differences in mean density and temperature. The respective sound speeds, wave constants, air 
densities, and Mach numbers are denoted by q and C 2 , fcj and kj, P\ and P 2 , and A/j and M 2 , 
where fcj=<a/q, £2 —col cj, Mi = Vj/q,and M2= Z V2^ C 2- 


In a hard-wall circular duct, the general solution to equation (2) is obtained as 

Here J m is the Bessel function of order m, p mn the nth zero of J' m (x), and A mn and B mn constant 
coefficients. The wave constants k^ n and k^ n correspond to the mode propagations, respectively, in 
the positive (to the right) and the negative (to the left) directions, and are given by 
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The integer n here is called the radial mode number, and the pair ( m,n ) is used to represent a single 
duct mode. 

Consider the incident wave of a single mode, say ( m ,£ ), 

< 5 ) 

V a ) 

This wave is incident from x=-°° and propagating toward duct termination, as shown in figure 1. 
Upon arriving at the duct termination, it will be partly reflected back into the duct and partly radiated 
out of the duct. In general, the reflected wave contains many radial modes, including propagating 
and attenuating modes, and is represented by 

rt = l " 1 ' 

Here R™ n is the conversion coefficient for the (m,^)mode incident and the ( m,n ) mode reflected. 
The reflection problem is completely solved by determining this coefficient for all values of n. The 
Wiener-Hopf technique yields radiation solutions in terms of the far field, which is represented by 

(7) 

* 2 K 

Here R is the radial distance from the center of the duct termination and 6 the polar angle measured 
from the x-axis (duct axis), as shown in figure 1. (R,d,<(>) are spherical coordinates. The complex 
factor f m ((0) is called the amplitude gain function, which provides the far field directivity of 
radiation. The phase A(k 2 ,M 2 ,R) of the far field depends on M 2 as well as on k 2 and R. The 
radiation problem is completely solved by determining f m e(0) and A(k 2 ,M 2 ,R)- 

This analysis with a single mode incident can be extended to accommodate incident waves 
composed of many modes in a straightforward manner. 

3. WIENER-HOPF FORMULATION 

The Wiener-Hopf technique involves extensive mathematical manipulation in the Fourier 
transform space. The Fourier transform of p(r, x) is given by 
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^(r,a) = -j==£^ p(r,x)c tax dx 


and p(r, x) is restored by the inverse transform 

I r 


p(r,x) = —j== I 0(r,a)e iajc da 


(9) 


In the process of the Wiener-Hopf formulation, various parameters and functions are defined and 
derived as follows: 
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In deriving K(a), we have used the condition of the continuity of the acoustic pressure and 
acoustic displacement at the mean flow mismatch at r = a, for x > 0. K(a) is factorized into 
two: one is analytic in the upper half plane (+) and the other in the lower half plane (-) as 
K(a)—K + (a) ■ K_(a). These factors will be included in the final solutions with arguments 
representing physical quantities. For the present problem, the factorization is obtained not in 
closed forms but in integral representations as follows: 


Lo g K +{y+) =±:\c + 

L° g K.^=sLj c _ 


Log K(a) ia 

(15) 

a-y+ 


Log *(«> da 

(16) 




Here the integral path is from -<» to +<» near the real axis, and the argument >’ + (y_ ) is located 
above (below) the respective integral path, as illustrated in figure 2. 
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Forgoing details of the formulation (ref. 5), the final results are presented here. The conversion 
coefficient for the reflection is obtained as 
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The symmetry between the radial mode numbers £ and n is salient, implying that the result satisfies 
the reciprocity principle (ref. 6), which can be used to infer the conversion coefficients of a 
nonpropagating mode. 

For the radiation, the phase is obtained as 
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Here the modified coordinates R' and 0' are defined as 


R'= (r 2 H — — — j , tan0'=^/l-M 2 2 tan0 

1 - M 0 Z 


( 21 ) 


5 


H$(x) is the Hankel function of the first kind, and H$ (x) is its derivative with respect 
to x. 

Comments are made on the expression of f m £(9) for two limiting cases. First, as can be 
shown readily, when 9 becomes zero, the quantity in the first square bracket will be infinite except 
for the case of m = 0. In other words, the radiated field is zero for 9 = 0 except for the radiation of 
axisymmetric modes. Second, as will be seen later, when its argument approaches -k+ n , K + in the 
second square bracket varies as - r/(0)] -1 . It follows that the radiated field will be zero for the 
angle satisfying 


Ti(9)=k+ n or cos0'=M 2 +(l-M 2 2 )^wL, for n±l 


( 22 ) 


However, for the angle corresponding to the incident propagation constant k^t , the radiated field 
is nonzero, because the term adjacent to K + becomes zero in this limit. In fact, the radiated field 
reaches the maximum in this limit. These findings are all familiar for cases of no mean flows. We 
will also see later that if its argument approaches —ky /(l+My), then K + varies in such a way as to 
compensate the term involving the square root in the same square bracket to maintain the amplitude 
gain finite. 


A remark should be made on the constant T ^ . As My and Af 2 both become zero, this constant 
becomes zero, but the expressions for the conversion coefficients and the radiation directivity remain 
correct, and finite when evaluated as a limit. 


4. NUMERICAL EVALUATION 


The integrals in equations (15) and (16) cannot be carried out analytically, and thus we will 
employ a semi-numerical method. To this end, all the variables are made dimensionless by multi- 
plying or dividing with the duct radius a. For notation simplification, the sub- or superscript m will 
often be dropped. 

K(a) satisfies all the conditions required for its factorization by the integration. Nevertheless, 
the integrand possesses singular points in the vicinity of the integral paths. These singular points 
arise as branch points, simple poles, and zeroes of K(a). It should be emphasized that there are no 
other singularities near the integral paths. The branch points are located at a= qj. We adopt a rule 
for determining phase around these branch points, as illustrated in figure 3. For example, consider 
(a-q[). Its phase is 180 degrees for its real part less than zero, and changes clockwise to zero as 
the real part becomes positive. On the other hand, the phase of (a-qy) is -180 degrees for the real 
part less than zero, and changes counterclockwise to zero as the real part becomes positive. This rule 
should be strictly observed for the integrations. 

The simple poles of K(a ) occur at zeroes of I' m (yj a), which is included through Wy (a) as in 
equations (1 1) and (14). Theses zeroes correspond to the wave constants of duct modes, and one can 
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show that the simple poles of K{(X) are located at a = v* s . Note that v„ (=-k nm ) is above 
the respective integral path, and v,7 (=-&7 ;/I ) is below the respective integral path, as shown in 
figure 2. K(a ) can also possess zeroes near the integral path if q+ > <72 • The zeroes are located 
between OC=q 2 and a=q\ , and above the integral path. The number of zeroes equals that of 
simple poles between a=q 2 and a=q\ , or can be less by one. The zeroes are denoted by z n , for 
n = 1, 2,... n 0 , with n 0 being the number of zeroes. These zeroes are ordered such that zj is the 
smallest and z no the largest. 

The imaginary parts of all the singular points are related to 4** ( k ). As the latter tends to zero, all 
the singular points approach the real axis, and the integral paths are then indented, as shown in 
figure 2(b). 

Consider the integral 



Log K(a) 

a-y 


da 


(23) 


This integral is divided, for convenience, as 

I = R_ + R++B_ + B + +S_+S + +Z+Y+N 

R’s are the contribution from the integration over larger arguments as 
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where x+ > j 

B’s are the contribution from the integration over small intervals containing the branch points: 
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J/h 


t >2 Log K{a) 
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B + = [*+ Lo 0 K(a) ^ * w ith b\ <q\ <l >2 (28) 

Jb x a-y 

S ' s are the contribution from the integration over small intervals containing the simple poles: 

«c 

S+ = ^ S~ , n c being the largest propagating radial mode number (29) 

n = l 
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Z’s are the contribution from the integration over small intervals containing the zeroes of K(a ): 


z=iz„ 

n=I 


_ rz „2 Log K(a) , 

Z n = I da , with z„i <z n < z n2 

J z n i a-y 


(32) 

(33) 


Y is the contribution from the interval containing the pole at a =y: 


K ± =f 

“ J y\ 


yi Log K(a) 


a-y 


da , with >1 < y < y 2 


(34) 


where the plus (negative) sign indicates that the pole is above (below) the integral path. This integral 
needs to be evaluated only if K(a) is free of singularity and zeroes within the integral limit. Other- 
wise, it should belong to one of B, S, and Z because there are no singular points other than those 
involved in B, S, or Z. 


Finally N is the integral over all the remaining intervals. There is no singularity at all in these 
intervals, and thus the integration can be carried out numerically. 


K{a ) approaches unity as |a| becomes large; that is, 
expanded for a large a , as 


Lim K(a) = 1, and thus K(a ) can be 
a— »±©o 


K{a) sl + ^ + 4 
a a 1 


(35) 


where the expansion coefficients A\ and A 2 can be readily obtained. With this substitution to 
equations (25) and (26), one obtains 


R± 
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where 


d=^-^- 


( 37 ) 


Error limits and the expansion coefficients are used to determine the integral limits x± ■ 
The integral near the branch points can be replaced by 

Log Q±(a) A „ 1 fLog(a-^f) 


J a-y 2 J a- 


■da 


where 


Q±{a) = K(a)-Ja^qf 

Q+(a ) are free of singularity and zeroes within the respective chosen integral limits. 
The simple poles of K(a) are separated as follows: 


? ± = f LogLj.(oQ i- Logta 

n J a-y a- 


Log (a-v~) 


da 


where 


-y 


L+(a)=K(a)-(a-v~) 


L+(a) are free of singularity and zeroes within the respective integral limits. 
The zeroes of K(a ) are similarly separated as 

_ r Log^(a) da + f Log(«-^) dn 

J a-y J a-y 


where 


U(a) = 


K(a) 

CC — Zn 


( 38 ) 


( 39 ) 


( 40 ) 


( 41 ) 


( 42 ) 


( 43 ) 


U(a) is free of singularity and zeroes within the respective integral limit. 
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Now consider the integral 


r*LogG(a) da 
J a a-y 


( 44 ) 


Here G(a) is free of singularity and zeroes within the integral limit, and thus represents Q±(a) in 
equation (38), L±(a) in equation (40), U(a ) in equation (42), or K(a) in equation (34) if K(a ) is 
free of singularity in that region. This integral can be written as 


H ( y) _ rfc Log[G(a)/G(y)] 
'a a-y 


da + Log G(y) 


rb da 
Ja a-y 


(45) 


As a approaches y, one readily obtains 

Log [G(a)/G(y)] _ G\y ) 
a-y G(y) 


(46) 


Thus, the first integral in equation (45) does not involve any singularity, and can be easily evaluated 
whether y is within the integral limit or not. The integral contained in the second term yields 


r o aa 

I = Log 

Ja a—y+ 

= Log 


(u „ > 

b-y± 

\ a -y±> 

' b-y ± ' 

y± ~ a 


, for y±<a or y±> b 
+ in , for a < y+< b 


(47) 


where (±) signs are used to indicate that the simple pole is above (+) or below (-) the integral path, 
which is indented around the pole like u for (+) and n for (-). 

Consider the integrals 


flr = f 0 Log(a-f-) da 
a-y_ 

(48) 

ctLo g (a-^) ia 
Ja a-y+ 

(49) 
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(51) 
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Here, the singular point at a= £ is within the integral limit, that is, a < £+ < b. The subscript on 
t, and y is used to indicate that the pole is above (+) or below (-) the integral path. The indented 
integral paths are shown in figure 4. The second integral in each of equations (38), (40), and (42) is 
identified with one of the integrals in equations (48) — (5 1 ). These integrals are evaluated as follows: 

For y± < a. 


f \ 2 

*2$ = +h rLog - ^+^-\{Log(b-y ± )) 2 +{Log(t; ± -y ± )) 2 

d-y± ) 3 2 L J 


- Log(£+ -;y+)- Log( a-y ± )+^ 
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^a-y ± V ' 
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Z±-y± 

b—y± 
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and for >>+ > b. 


Q+ = +/7rLog 
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+ Log(y+ - £+ ) ■ Log(>'+ -£)-£- 2 
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y±-b 
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y±z3± 
y y±~ a J 


where the upper (lower) sign of in corresponds to (£_). 
For a <y±<b. 


(52) 


(53) 


Qz = ± in Log()'_ -a) + in [- Log(;y_ - £+ ) + Log(j_ -<£+)] 

- ^(Log(Z? - £+ )) 2 + (Log(y_ - a)) 2 + Log(b-£+)-Log(b-y_) 
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2 ;= 1 j 2 


r y--J±) J J y--Z± 


b~£± 
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y--a 


Q+ = ± /7T Log(>' + -a) + i^[ Log(>> + - £+) + Log(y+ - £+)] 

-2 (Log(Z?-<^ ± )) 2 +(Log(y + - a)) 2 ] + Log(fc -£ ± ) • Log(6 - y+) 
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(54) 


(55) 
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All the integrals involving singular points are analytically evaluated according to equations (47) 
and (52)— (55). For large arguments extending to the infinite, the integral is evaluated according to 
equation (36). The rest involves finite integrals with well-behaved integrands, and thus can be 
numerically evaluated without any difficulty. 

A remark will be made on the results in equations (54) and (55), particularly on the second term 
on the right. This term vanishes if y and£ are on the same side (above or below) with respect to 
the integral path. However, if y and£ are separated by the integral path, the second term equals 
±2m'Log(y+ -<*+) and diverges as y approaches £. Of the -factors contained in the results in 
equations (17) and (19), only K + (-T](0)) can have y and £ separated by the integral path. Its 
argument is above the integral path according to equation (15). As 6 varies, the argument can 
approach a singular point located below the integral path. For this case, the second term above is 
given by 2mLog(-77(0)-£_), where is v~ or q[ . Inspecting equations (15), (23), (24), (38), 
and (55), one obtains, for -rj(6) close to , 

K + (-tK0))oc-j L_ (56) 

7(0) -<?i 

Similarly, from equations (15), (23), (24), (40), and (55), one obtains, for -rj(6) close to v n , 

M-m0))~(-m0)-v;) _1 ( 57 ) 

It follows that, as rj(8) approaches k ^ n , then given in equation (19) tends to be zero for 

n^l and will reach the maximum for n-i. Also, one can see that f m g(9) will remain finite as rj(6) 
approaches k\/( 1 + M \) , as discussed earlier. 


CONCLUDING REMARKS 


An analysis for evaluation of the Wiener-Hopf solution was presented for sound radiation from 
an unflanged circular duct with mean flows. This analysis was initially developed for radiation of 
spinning modes in conjunction with aircraft inlet noise control studies. We have a well-working 
computer code available for such radiation. However, while generating numerical results for the 
benchmark problems, we learned that the code needed to be refined for radiation of the 
axisymmetric modes. 
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Figure 1. Sound radiation from unflanged circular duct with flow. 
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(b) 

Figure 2. (a) Integral paths C+ and C_ and singular points ( V 3 corresponds to an attenuating 
mode); (b) indentation of the integral paths as fa (k) — > approaches 0 . 
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Figure 4. Integral path indentation for integrations involving two singular points at a = y and c, 
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